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ABSTRACT 

Point-to-point correlation functions of hadron currents in the QCD vacuum are calculated 
on a lattice and analyzed using dispersion relations, providing physical information down to 
small spatial separations. Qualitative agreement with phenomenological results is obtained in 
channels for which experimental data are available, and these correlation functions are shown 
to be useful in exploring approximations based on sum rules and interacting instantons. 
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One of the major challenges of lattice QCD is to provide as much insight as possible 
and to address as wide a range of experimental observables as possible using the Euclidean 
observables which are amenable to lattice calculations. In the past, the primary emphasis has 
been on correlation functions at large Euclidean separations. For example, the asymptotic 
decay of correlation functions between widely separated hadronic currents is used to measure 
hadron masses,^ distantly separated hadronic sources and sinks are used to filter out hadronic 
ground states, from which ground state observables, wave functions, and form factors are 
obtained,^ and spatially well-separated hadronic sources are used at finite temperature to 
extract analogous screening masses and wave functions.^ In these large distance calculations, it 
has been useful to integrate one of the currents over space to project onto a specific momentum 
and non-local sources are usually employed to increase the overlap with the desired states. 

This present work addresses complementary QCD physics at short and intermediate dis- 
tances, as well as large separations, by calculating the point-to-point Euclidean correlation 
functions R{x) = (0 \TJ{x)J{0) \ 0) where J is one of the point hadron currents listed in Ta- 
ble I. Physically, these Euclidean correlation functions correspond to space-like-separated or 
equal-time correlation functions describing the virtual propagation of quarks or hadrons. They 
complement bound state hadron properties in the same way that nucleon-nucleon scattering 
phase shifts provide detailed information about the spatial distribution of spin-spin, spin-orbit 
and tensor forces complementary to that provided by properties of the deuteron. These corre- 
lation functions have been determined phenomenologically in some channels by using disper- 
sion relations to analyze hadron production and r-decay experiments,^ and have been studied 
extensively using QCD sum rules and the interacting instanton approximation.^"^ Hence, in 
this work, we report what we believe to be the first exploratory lattice QCD calculation of 

2 



the behavior of point-to-point meson and baryon correlation functions, from short distances 
where the physics is perturbative to large distances where it is highly non-perturbative. We 
demonstrate promising agreement with experimentally measured results where available, show 
how to analyze the lattice measurements using dispersion relations, and indicate the potential 
more extensive calculations have to explore, test, and refine approximations. 

The lattice calculations were performed on a 16^ x 24 lattice in the quenched approxi- 
mation at inverse coupling ^jg^ = 5.7, corresponding to a physical lattice spacing defined by 
the proton mass of approximately 0.168 fm. This inverse coupling constant is large enough 
to give a semi-quantitative approximation and has the significant advantage that light quark 
propagators are available for point sources. Hence, propagators with point sources calculated 
by Soni et al.^ for 16 configurations with five values of the quark mass were used to calculate 
the meson and baryon correlation functions listed in Table I. For the benefit of non-specialists, 
we convert the hopping parameters n = (0.154, 0.160, 0.164, 0.166, and 0.168) and critical 
value Kc = 0.1692 to the intuitively useful bare quark mass m = (1/2k - 1/2kc) a-^ which 
takes on the values (351, 199, 110, 67, and 25) MeV. Extrapolation to the value of k, at which 
the pion acquires its correct mass corresponds to extrapolation to m = 8 MeV. In order to 
calculate the ratio of interacting to free correlation functions, corresponding free lattice prop- 
agators at a very small quark mass, ma = 0.05 were calculated exactly on a (48)^ lattice, 
where the lattice volume was chosen large enough to eliminate finite volume effects for spatial 
separation less than 2 fm and the value of ma was chosen small enough to produce deviations 
from the massless result less than one percent at 1 fm. Because the propagators had hard- 
wall boundary conditions in the time direction, all correlation functions were calculated on 
the central time slice containing the source. 
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Several lattice artifacts had to be eliminated to obtain a physical approximation to the 
continuum correlation functions. Significant anisotropy is introduced in the rotationally in- 
variant continuum correlation functions by the Cartesian lattice. Although at very short 
physical distances the anisotropics in the numerator and denominator cancel by asymptotic 
freedom and at a sufficiently large number of lattice spacings the granularity of the lattice 
becomes negligible, for the present lattice spacing there remains an intermediate region in 
which some correlation function ratios display substantial anisotropy. As discussed and jus- 
tified in detail in a subsequent article,^ since the lattice points far away from the Cartesian 
directions agree with the continuum result in the non-interacting case and should be the most 
reliable in the interacting theory, the lattice correlation functions R{x) in this work are taken 
from sites x within an angular cone surrounding the diagonal directions d = (n, n, n) such 
that d ■ X > 0.9. Also, although in principle one would like to normalize the ratios at an 
infinitesimally small separation, we have had to normalize our results at the first non-zero 
diagonal separation (1, 1, 1), corresponding to a physical separation of ~ 0.29 fm. 

A second lattice artifact arises from the contributions of periodically repeated images of 
the physical sources on the finite lattice having periodic boundary conditions. As discussed in 
detail in Ref. [10], it is straightforward to perform a self-consistent subtraction of the image 
contributions, and we have applied and verified this image correction to all the correlation 
functions presented in this work. 

Finally, because it is impractical to calculate quark propagators at a quark mass light 
enough to produce a physical pion, it is necessary to perform a sequence of calculations at a 
series of heavier quark masses and extrapolate to the pion limit. Although we do not know 
the proper functional form for the extrapolation in the chiral limit, the extrapolation is in fact 
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innocuous because we only have to extrapolate data at ruq = 351, 199, 110, 67, and 25 MeV 
down to 8 MeV. A worst case example will be shown below in Fig. 1. Thus, we believe the 
technical aspects of relating the present lattice calculations to the physical continuum theory 
are under reasonable control, and that our results provide a meaningful first comparison of 
quenched QCD with phenomenological and model results. 

The correlation functions we considered are listed in Table I. For convenience, we have 
contracted the vector indices in the vector and axial channels. In addition, we inserted a 
factor of Xfj^ji^ in the baryonic correlation functions before taking the Dirac trace in order 
to project out the component which is stable in the massless quark limit. For a local field 
theory such as QCD, a two-point function is uniquely characterized by its absorptive part in 
momentum space (up to a polynomial) through the dispersion relation 



where the upper, middle and lower terms in brackets refer to scalar and pseudoscalar mesons, 
vector and axial vector mesons and baryon channels respectively and terms which vanish 
for the correlation functions in Table I have been omitted. Phenomenologically, we expect 
f{s) has two major contributions, a resonance piece and a continuum piece, parameterized 
as f{s) = X'^5{s — M^) + fp{s)9{s — sq), where M is the resonance mass, A is the coupling of 
the current to the resonance state and fp{s) is the perturbative contribution. Since QCD is 
asymptotically free, the /p(s)'s are well-approximated by the corresponding free correlation 
functions with massless quarks, also listed in Table I. This type of parameterization is widely 
used in QCD sum rule calculations. 




(1) 
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An inverse Fourier transform of Eq. (1) defines the following phenomenological correlation 
functions in coordinate space as a function of M, A and sq, where the polynomial P{q^) only 
contributes at the point x = and can be ignored for finite x: 

3M^x-^ I ^/ 1 \ J 1 1 ^ (2) 



Here, m refers to the upper terms in brackets for scalar and pseudoscalar mesons and the 
middle terms in brackets for vector and axial vector mesons, and b refers to the lower terms 
in brackets for baryon channels. In practice, we always normalize the phenomenological R{x) 
by the corresponding massless results R'o'ix) oc x' and Rl{x) (XX ^. Note that although 
we have not included lattice renormalization constants for non-conserved currents, they will 
not affect the shape of the normalized ratio of correlation functions, but only the value of the 
fitted parameters. Also, since the results below are sensitive to the presence of the continuum 
term in the fitting, it is clear that the lattice calculation is simultaneously providing short 
and long distance physics information. 

In Fig. 1 we show the complete lattice data and the three-parameter fit in the pseudoscalar 
meson (pion) channel for each of five values of the bare quark mass. The result at the physical 
pion mass is obtained by binning the lattice data in two-lattice- unit bins and extrapolating the 
binned data with the results shown by the solid circles. The striking result in this channel is the 
extremely rapid rise in the correlation function ratio, necessitating a log plot, which arises from 
the strong attraction and corresponding light pion mass. The lattice result agrees qualitatively 
with Shuryak's phenomenological estimate,^ denoted by the dot-dashed line, based on the 
value = (480 MeV)^ and the fact that the peak is proportional to X-^/rri^ ~ fn/iT^q explains 
the particularly large quark mass dependence in this channel. Detailed treatment of the 
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extrapolation, error analysis of fitted parameters, and lattice renormalization corrections will 
be deferred to the longer paper, ^ and we only show extrapolated results and phenomenological 
fits for all other channels, where the quark mass dependence is much weaker. 

Figure 2a displays the vector meson (rho) channel result. As emphasized by Shuryak,^ 
the salient feature in this channel is the fact that although the free correlator falls four orders 
of magnitude between 0.3 and 1.5 fm, the ratio is nearly one over the whole range, and our 
lattice result is consistent with his phenomenological analysis of e+e" — > even number tt's, 
denoted by the dashed curve. The result in the axial meson (Ai) channel, shown in Fig. 2b, 
is qualitatively similar to the phenomenological analysis of r — > Stt decay^ denoted by the 
dashed line, although finite lattice effects render it difficult to reproduce the rising tail due to 
mixing with the pion. The result for the scalar meson channel, for which the extrapolation 
was more problematic than any other, is shown in Fig. 2c. Note that in both the axial and 
scalar channels, Eq. (2) did not produce reasonable physical parameters when fit to the data, 
so the solid curves are smooth curves to guide the eye in these cases. For comparison, the 
predictions of the interacting instanton approximation^ for mesons using a Pauli-Villars cutoff 
Apv = 130 MeV are shown by dotted lines, and in each case are in qualitative agreement with 
the lattice results. Comparison of calculations on the same lattice with "cooled" configurations 
which are in progress will be particularly instructive in this connection. 

Since there are no phenomenological data in the baryon channels, in Fig. 3 we have 
compared nucleon and delta correlators with the sum rule calculations of Belyaev and loffe^ 
(dashed lines) and Farrar et alJ (dot-dashed lines). Given the agreement with phenomenol- 
ogy in the meson channels and drastic differences in sum rule results for the delta, lattice 
correlation functions can play a useful role in testing and refining sum rule approximations. 
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In summary, we believe these exploratory calculations demonstrate the feasibility and utility 
of calculating vacuum correlation functions on a lattice and motivate definitive calculation on 
larger lattices with 6/(7^ > 6. 

It is a pleasure to thank Edward Shuryak for stimulating discussions at the Aspen Center 
for Physics and elsewhere and for making his data available to us prior to publication. We 
also thank Amarjit Soni for making his point propagators available to us, and the National 
Energy Research Supercomputer Center for Providing Cray-2 computer resources. 
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FIGURE CAPTIONS 

Fig. 1: Ratio of interacting to free correlation functions measured at five quark masses (open cir- 
cles), binned data extrapolated to the physical pion mass (closed circles), three-parameter 
fits (solid lines) and the phenomenological results of Ref. [1] (dot-dashed line). 



Fig. 2: Extrapolated ratio of meson correlation functions (closed circles) and fits (solid curves) 
as in Fig. 1. Dashed lines denote phenomenological results^ and dotted lines show the 
interacting instanton approximation.^ 



Fig. 3: Extrapolated ratio of baryon correlation functions (closed circles) and fits (solid curves) 
as in Fig. 1. Dashed and dot-dashed lines denote sum rule results from Ref. [6] and [7], 
respectively. 
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Table I: Hadron currents and correlation functions. 


Channel 


Current 


Correlator 




Vector 


,/,, = Ti,'y,,d, 


(O \TJ„(x)Jn(0)\ O) 


1 

127r2 


Axial 






1 

127r2 


PspiidoscalaT 




lTJPfx)J^'fO)l 0\ 

\ V7 1 -1. I %AJ ) tJ \ / 1 / 


3s 
87r2 


Scalar 


J* = ltd 


(O TJ*(a;)P(0) O) 


3s 
87r2 


Nucleon 




^Tr ((0|TJ^(a;)J^(0)|0>x^7^) 


S2 
6477^ 


Delta 


= eabc [«"C7^w'] 


((0 |TJ^(x)J^(0)|0>x.7l 


3^2 
2567r4 
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